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We present the results of a muon-spin relaxation (p^SR) investigation into magnetic ordering 
in several families of layered quasi-two-dimensional molecular antiferromagnets based on transition 
metal ions such as 5* = | Cu^^ bridged with organic ligands such as pyrazine. In many of these ma- 
terials magnetic ordering is difficult to detect with conventional magnetic probes. In contrast, p^SR 
allows us to identify ordering temperatures and study the critical behavior close to Tn . Combining 
this with measurements of in-plane magnetic exchange J and predictions from quantum Monte Carlo 
simulations we may assess the degree of isolation of the 2D layers through estimates of the effective 
inter- layer exchange coupling and in- layer correlation lengths at Tn. We also identify the likely 
metal-ion moment sizes and muon stopping sites in these materials, based on probabilistic analysis 
of the magnetic structures and of muon-fluorine dipole-dipole coupling in fluorinated materials. 

PACS numbers: 76.75.+i, 75.50.Xx, 75.10.Jm, 75.50.Ee 



I. INTRODUCTION 

The S = ^ two-dimensional square-lattice quantum 
Heisenberg antiferromagnet (2DSLQHA) continues to 
be one of the most important theoretical models in 
condensed matter physics^. Experimental realizations 
of the 2DSLQHA in crystals also contain an interac- 
tion between planes, so that the relevant model de- 
scribing the coupling of electronic spins Si gives rise 
to the Hamiltonian^ 

where J {J±) is the strength of the in- (inter-) plane 
coupling and the first (second) summation is over 
neighbors parallel (perpendicular) to the 2D x?/-plane. 
Any 2D model {J± = 0) with continuous symmetry 
will not show long-range magnetic order (LRO) for 
T > due to a divergence of infrared fluctuations^'^. 
However, layered systems approximating 2D models 
{J_\_ 7^ 0) will inevitably enjoy some degree of inter- 
layer coupling and this will lead to magnetic order, 
albeit at a reduced temperature due to the influence 
of quantum fluctuations. Quantum fluctuations are 
also predicted to reduce the value of the magnetic mo- 
ment in the ground state of the 2DSLQHA to around 
60% of its classical value^, and this reduction is often 
seen in the ordered moments of real materials. In lay- 
ered materials that approximate the 2DSLQHA, the 
measurement of the antiferromagnetic ordering tem- 
perature Tn is often problematic due not only to this 
reduction of the magnetic moment, but also to short- 
range correlations that build up in the quasi-2D lay- 
ers above Tn- These correlations lead to a reduction 



in the size of the entropy change that accompanies 
the phase transition, reducing the size of the anomaly 
in the measured specific heat^. We have shown in a 
number of previous cases that muon-spin relaxation 
(p+SR) measurements do not suffer from these effects 
and therefore represent an effective method for detect- 
ing magnetic order in complex anisotropic systems^^ . 

The rich chemistry of molecular materials allows 
for the design and synthesis of a wide variety of 
highly-tunable magnetic model systems^^. Magnetic 
centers, exchange paths and the surrounding molec- 
ular groups can all be systematically modified, al- 
lowing investigation of their effects on magnetic be- 
havior. In particular, the existence of different ex- 
change paths along different spatial directions can re- 
sult in quasi-low-dimensional magnetic behavior (i.e. 
systems with magnetic interactions constrained to act 
in a two-dimensional plane or along a one-dimensional 
chain). Such systems have the potential to better ap- 
proximate low-dimensional models than many tradi- 
tional inorganic materials. In addition, these molecu- 
lar materials can have exchange energy scales of order 
J//cb ~ 10 K which are accessible with typical labo- 
ratory magnetic fields^ of 5 ~ 10 T allowing an addi- 
tional avenue for their experimental study. This con- 
trasts with typical inorganic low-dimensional systems 
where the exchange is found to be J//cb ~ 1000 K and 
fields of B ^ 1000 T would be needed to significantly 
perturb the spin system. 

It has also been shown^^'^^ that a small XY-\ike 
anisotropy exists in some molecular materials. Al- 
though J± is the decisive energy scale for the mag- 
netic ordering, this anisotropy has been shown to have 
an influence on the ordering temperature^^ and deter- 
mines the shape of the low-field B-T phase diagram 
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Figure 1. Bridging ligands used in the compounds de- 
scribed in this paper: (a) pyrazine (N2C4H4, abbreviated 
pyz); and (b) (ii) pyridine- AT-oxide (C5H5NO, abbreviated 
pyo). Dashed hues indicate where the hgands bond to 
other parts of molecular structures. 



of these systems-^^. 

Several classes of molecular magnetic material 
closely approximate the 2DSLQHA model and in 
this paper we report the results of p+SR measure- 
ments performed on several such materials. These 
systems are self- assembled coordination polymers, 
based around paramagnetic ions such as Cu^+, linked 
by neutral bridging ligands and coordinating anion 
molecules. Our materials are based on combinations 
of three different ligands: (i) pyrazine (N2C4H4, ab- 
breviated pyz) and (ii) pyridine- TV-oxide (C5H5NO, 
abbreviated pyo), both of which are planar rings; and 
(iii) the linear bifluoride ion [(HF2)~], which is bound 
by strong hydrogen bonds F - R - F. The pyz and pyo 
ligands are shown in Fig. 1. Specifically, we investi- 
gate the molecular system [M(HF2)(pyz)2]X, where 
= Cu^+ is the transition metal cation and X~ 
is one of various anions (e.g. BF4 , CIO4 , PFg etc.). 
We also report the results of our measurements on 
other quasi-2D systems. First, [Cu(pyz)2(pyo)2] 1^2, 
with Y~ = BF4 or PF^ , in which pyo ligands bridge 
Cu(pyz)2 planes. Then, the quasi-2D non-polymeric 
compounds [Cu(pyo)6]^2, where Z~=BF^^ CIO3 or 
PF^ are examined. We also investigate materials in 
which either Ni^^ {S = 1) or Ag2+ {S = ^) form the 
magnetic species in the quasi-2D planes rather than 
Cu2+. 



This paper is structured as follows. In Sec. II we 
outline the p+SR technique and describe our exper- 
imental methods. The [Cu(IIF2)(pyz)2]X family of 
materials is then discussed in Sec. Ill, where muon 
data are used to determine Tn and critical parame- 
ters. Muon-fluorine dipole-dipole oscillations in the 
paramagnetic regime are found for these materials 
which we use, in conjunction with dipole field sim- 
ulations, to investigate possible muon sites and con- 



strain the copper moment. In Sec. IV we explore 
the related 2D system [Cu(pyz)2(pyo)2] y^2- Sec. V 
details measurements of [Cu(pyo)6]^2- Sec. VI ex- 
amines a highly two-dimensional silver-based molec- 
ular material, Ag(pyz)2(S2 08). Finally, data from 
the [Ni(HF2)(pyz)2]X (X" = PFg , SbF^ ) family of 
molecular magnets is presented in Sec. VIL 



II. EXPERIMENTAL DETAILS 

Zero- field (ZF) p+SR measurements were made on 
powder samples of the materials at the ISIS facil- 
ity, Rutherford Appleton Laboratory, UK using the 
MuSR and EMU instruments and the Swiss Muon 
Source (SpS), Paul Scherrer Institut, Switzerland us- 
ing the General-Purpose Surface-Muon (GPS) instru- 
ment and Low- Temperature Facility (LTF). For mea- 
surements at temperatures T > 1.8 K powder samples 
were packed in a 25 pm Ag foil packet and mounted on 
a Ag backing plate. For measurements at T < 1.8 K 
the samples were mounted directly on an Ag plate and 
covered with a 12.5 pm Ag foil mask. Ag is used since 
it has only a small nuclear magnetic moment, and so 
minimizes the background depolarization of the muon 
spin ensemble. 

In a p+SR experiments^, spin-polarized positive 
muons are stopped in a target sample. The posi- 
tive muons are attracted to areas of negative charge 
density and often stop at interstitial positions. The 
observed property of the experiment is the time evo- 
lution of the muon-spin polarization, the behaviour 
of which depends on the local magnetic field at the 
muon site. Each muon decays with an average life- 
time of 2.2 ps into two neutrinos and a positron, the 
latter particle being emitted preferentially along the 
instantaneous direction of the muon spin. Recording 
the time dependence of the positron emission direc- 
tions therefore allows the determination of the spin 
polarization of the ensemble of muons. In our experi- 
ments, positrons are detected by detectors placed for- 
ward (F) and backward (B) of the initial muon polar- 
ization direction. Histograms N-p{t) and NB{t) record 
the number of positrons detected in the two detectors 
as a function of time following the muon implanta- 
tion. The quantity of interest is the decay positron 
asymmetry function, defined as 



A{t) = 



Np{t) - aNsit) 
NF{t) + aNB{t) 



(2) 



where a is an experimental calibration constant. The 
asymmetry, A(t), is proportional to the spin polariza- 
tion of the muon ensemble. 

A muon spin will precess around the local mag- 
netic fields at its stopping site at a frequency u = 
7p5/27r, where the muon gyromagnetic ratio 7^ = 
27r X 135.5 MIIzT"^ In the presence of LRO in a 
material we often measure oscillations in A{t). These 
result from a significant number of muons stopping 
at sites with a similar internal field, giving rise to a 
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Figure 2. Example A{t) spectra with fits above (T = 25 K) 
and below (T = 0.35 K) the magnetic transition temper- 
ature Tn = 1.4 K for [Cu(HF2)(pyz)2]BF4. Note the ap- 
proximately equal initial asymmetry, and the oscillations 
and '|-tair observed in the ordered phase. The slow oscil- 
lation observed for T > Tn is due to muon-fluorine dipole- 
dipole oscillations. 

coherent precession of the ensemble of muon spins. 
Since the spins precess in local magnetic fields directed 
perpendicularly to the spin polarization direction, we 
would expect that, for a powder sample with a static 
magnetic field distribution, | of the total spin com- 
ponents should precess and the remaining ^ should 
be non-relaxing. The non-relaxing third of muon- 
spin components give rise to the so-called |-tail in 
A{t)^ whose presence therefore provides additional ev- 
idence for a static field distribution in powder sample. 
Taken together, these effects provide an unambigu- 
ous method for sensitively identifying a transition to 
LRO. An example of typical spectra above and be- 
low the magnetic ordering temperature is shown for 
[Cu(HF2)(pyz)2]BF4 in Fig. 2; the oscillations above 
Tn are characteristic of a quantum-entangled Fp state 
(see Sec. HID). 

III. [Cu(HF2)(pyz)2]X 

The synthesis of the [M(HF2)(pyz)2]X system^' ^^^'^^ 
represented the first example of the use of a bifluoride 
building block to make a three-dimensional coordi- 
nation polymer. This class of materials possesses a 
highly stable structure due to the exceptional strength 
of the bifluoride hydrogen bonds. The structure of the 
[M(HF2)(pyz)2]X system^^'-^^ comprises infinite 2D 
[M(pyz)2]^^ sheets which lie in the ab plane, with bi- 
fluoride ions (HF2)~ above and below the metal ions, 
acting as bridges between the planes to form a pseu- 
docubic network. The X~ anions occupy the body- 
center positions within each cubic pore. An example 
structure, for [Cu(HF2)(pyz)2]PF6, is shown in Fig. 
3. Samples are produced in polycrystalline form via 
aqueous chemical reactions between MX2 salts and 
stoichiometric amounts of ligands. Preparation de- 
tails for the compounds are reported in Refs. 14-16. 

In this section we consider those materials where 




Figure 3. The structure of [Cu(HF2)(pyz)2]PF6, as an ex- 
ample of the [M(HF2)(pyz)2]X series. Copper ions are 
joined in a 2D square lattice by pyrazine ligands to form 
Cu(pyz)2^ sheets; the 2D layers are joined in the third di- 
mension by HF^ groups, making a pseudocubic 3D struc- 
ture; and this structure is stabilised by a PF^ anion at 
the center of each cubic pore. For clarity, hydrogen atoms 
attached to pyrazine rings have been omitted, and only 
one PFr anion is shown. 



the M cations are Cu^+ 3d^ S = ^ centers. It is 
thought that the magnetic behavior of these mate- 
rial results from the 3dx2_y2 orbital of the Cu at the 
center of each CUN4F2 octahedron lying in the CUN4 
plane so that the spin exchange interactions between 
neighboring Cu^+ ions occur through the 5-bonded 
pyz ligands^^. The interplane exchange through the 
IIF2 bridges connecting two Cu^+ ions should be very 
weak as these bridges lie on the 4-fold rotational axis 
of the Cu 3dx2_y2 magnetic orbital, resulting in lim- 
ited overlap with the fluorine Pz orbit als. Therefore 
to a first approximation, the magnetic properties of 
[M(IIF2)(pyz)2]X can be described in terms of a 2D 
square lattice. 

Measurements for X~ = BF4 , CIO4 and SbFg 
were made using the MuSR spectrometer at ISIS, 
whilst PFe, AsFe, NbFe and TaFe were measured us- 
ing GPS at PSI. 



A. Long-range magnetic order 

The main result of our measurements on these sys- 
tems is that below a critical temperature Tn, oscilla- 
tions in the asymmetry spectra A{t) are observed at 
two distinct frequencies, for all materials in the series. 
This shows unambiguously that each of these materi- 
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Figure 4. Data and the results of fits to Eq. (3) for [Cu(HF2)(pyz)2]X magnets with tetrahedral anions X~ . From left 
to right: (a) and (d) show sample asymmetry spectra A{t) for T < Tn along with a fit to Eq. (3); (b) and (e) show 
frequencies as a function of temperature [no data points are shown for the second line because this frequency 1^2 was held 
in fixed proportion to the first, ui (see text)]; and (c) and (f) show relaxation rates as a function of temperature. In 
the iy{T) plot, error bars are included on the points but in most cases they are smaller than the marker being used. The 
solid line representing in (b) and (c) corresponds to the filled circles in the third column of graphs [(c) and (f)] for 
that component's relaxation, Ai, whilst the dashed line and unfilled circles correspond to 1^2 and A2, respectively. The 
filled triangles correspond to the fast relaxation A3. 



als undergoes a transition to a state of LRO. Example 
asymmetry spectra are shown in the left-hand column 
of Fig. 4 and Fig. 5. They were found to be best fitted 
with a relaxation function 

A{t) =Ao [pie~^^* cos{2Tivit + ^1) 

+P2e"^^^ cos(27rP2^it + (^^2) + Pse"^^'] 

+Abge-^^^^ (3) 

where Aq represents the contribution from those 
muons which stop inside the sample and Abg accounts 
for a relaxing background signal due to those muons 
that stop in the silver sample holder or cryostat tails, 
or with their spin parallel to the local field. Of those 
muons which stop in the sample, pi indicates the 
weighting of the component in an oscillating state with 
frequency z^i; p2 is the weighting of a lower- frequency 
oscillating state with frequency and p^, represents 
the weighting of a component with a large relaxation 
rate A3. All parameters were initially left free to vary. 
The second frequency was found to vary with temper- 
ature in fixed proportion to vi via V2 = ^2^1 for each 
material. The parameter P2 was identified by fitting 
the lowest-temperature A(t) spectra where Eq. (3) 
would be expected to most accurately describe the 
data, and subsequently held fixed during the fitting 
procedure. Phase factors (pi were also found to be 
necessary in some cases to obtain a reliable fit. The 
parameters resulting from these fits are listed in Ta- 
ble I, and data with fits are shown in Figs. 4 and 5. 
We also note here that the discontinuous nature of 



the change in all fitted parameters and the form of 
the spectra at Tn strongly suggest that these materi- 
als are magnetically ordered throughout their bulk. 

The frequencies and relaxation rates as a function 
of temperature extracted from these fits are shown in 
the central column of Figs. 4 and 5. The muon preces- 
sion frequency, which is proportional to the internal 
field in the material, can be considered an effective 
order parameter for the system. Consequently, fitting 
extracted frequencies as a function of temperature to 
the phenomenological function 



v{T) = uiO) 



T 



(4) 



allows an estimate of the critical temperature and the 
exponent /3 to be extracted. Our results fit well with a 
previous observation^ that the compounds divide nat- 
urally into two classes: those with tetrahedral anions 
X~ = BF4 , CIO4 and those with octahedral anions 
X~ = AFq . The tetrahedral compounds have lower 
transition temperatures Tn ^ 2 K, as compared to the 
octahedral compounds' Tn ^ 4 K; and the tetrahedral 
compounds also display slightly lower oscillation fre- 
quencies than their octahedral counterparts'". 

This difference has been explained in terms of 
differences in the crystal structure between the two 
sets of compounds. Firstly, the octahedral anions are 
larger than their tetrahedral counterparts. Secondly, 
the pyrazine rings are tilted by differing amounts with 
respect to the normal to the 2D layers: those in the 
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Figure 5. Example data and fits for [Cu(HF2)(pyz)2]X magnets with octahedral anions X~ . From left to right: (a), 
(d), (g), (j) and (m) show sample asymmetry spectra A(t) for T < Tn along with a fit to Eq. (3); (b), (e), (h), (k) and 
(n) show frequencies as a function of temperature [no data points are shown for the second line because this frequency 
1^2 was held in fixed proportion to the first, ui (see text)]; and (c), (f), (i), (1) and (o) show relaxation rates as a 
function of temperature. In the iy{T) plot, error bars are included on the points but in most cases they are smaller than 
the marker being used. The solid line representing ui in (b), (e), (h), (k) and (n) corresponds to the filled circles in the 
third column of graphs [(c), (f), (i), (1) and (o)] for that component's relaxation, Ai, whilst the dashed line and unfilled 
circles correspond to 1^2 and A2, respectively. 



octahedral compounds are significantly more upright. 
Since the Cu 3dx2-y2 orbitals point along the pyrazine 
directions, these tilting angles might be expected, 
to first order, to make little difference to nearest- 



neighbor exchange because such rotation is about a 
symmetry axis as viewed from the copper site. How- 
ever, it may be that the different direction of the de- 
localized orbitals above and below the rings through 
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Table I. Fitted parameters for molecular magnets in the [Cu(HF2)(pyz)2]X family. The first parameters shown relate 
to fits to Eq. (3), which allow us to derive frequencies at T = 0, Ui; probabilities of stopping in the various classes of 
stopping site, pi, in percent; and phases associated with fitting the oscillating components 0^. Then, the temperature 
dependence of z/^ is fitted with Eq. (4), extracting values for the Neel temperature, Tn, critical exponent [3 and parameter 
a. Finally, the quoted J/fe is obtained from pulsed-field experiments^, and the ratio of inter- to in-plane coupling, J±/J, 
is obtained by combining Tn and J with formulae extracted from quantum Monte Carlo simulations (see Sec. IIIB, and 
Ref. 7). Dashes in the A3 column for the SbFe and TaFe compounds indicate that no fast-relaxing component was used 
to fit those data. Dashes in the J/fe and J±/J columns for NbFe and TaFe indicates a lack of pulsed-field data for 
these materials. 



which exchange probably occurs, possibly in conjunc- 
tion with hybridization with the anion orbit als, result 
in an altered next-nearest neighbor or higher-order in- 
teractions, changing the transition temperature. 

Within the tetrahedral compounds, the difference in 
the weighting of the oscillatory component {pi -\- P2) 
in X~ = BF4 and CIO4 probably results from the 
difficulty in fitting the fast-relaxing component. Even 
with little change in the size of the oscillations, any 
error assigning the magnitude of this component will 
affect the proportion of the A{t) signal attributed to 
them. This difficulty is partly due to the resolution- 
limited nature of ISIS arising from the pulsed beam 
structure. In the octahedral compounds, we found 
that X~ = SbFg and TaFg did not have a resolvable 
fast-relaxing component, and consequently ps was set 
to zero during the fitting procedure. This is reflected 
by dashes in the ps and A3 columns in Table I. 

The fact that two oscillatory frequencies are ob- 
served points to the existence of at least two mag- 
netically distinct classes of muon site. In general we 
find that pi ^ p2 for these materials, making the prob- 
ability of occupying the sites giving rise to magnetic 
precession approximately equal. The weightings pi^2 
were found to be significantly less than the weighting 
Ps relating to the fast-relaxing site. This, in com- 
bination with the magnitude of the fast relaxation 
A3(T = 0) > 10 MHz, suggests that this term should 
not be identified with the ^-tail which results from 
muons with spins parallel to their local field. (If that 
were the case then we would expect {pi -\-p2)/p3 = 2, 
which we do not observe.) It is likely that each of 
the components, pi, p2 and pa, therefore reflect the 
occurence of a separate class of muon site in this sys- 
tem. We investigate the possible positions of these 
three classes of site in Sec. HIE. 

The temperature evolution of the relaxation rates 
Xi is shown in the right-hand columns of Figs. 4 
and 5. In the fast-fluctuation limit, the relaxation 
rates are expected^^ to vary as A oc A^r, where 

A = J^'^{{B — BqY) is the second moment of the 



local magnetic field distribution (whose mean is ^o) 
in frequency units, and r is the correlation time. In 
all measured materials, the relaxation rate Ai, cor- 
responding to the higher oscillation frequency, starts 
at a small value at low temperature and increases as 
Tn is approached from below. This is the expected 
temperature-dependent behavior and most likely re- 
flects a contribution from critical slowing down of fluc- 
tuations near Tn (described e.g. in Ref. 18). In con- 
trast, the relaxation rate A2 (associated with the lower 
frequency) starts with a higher magnitude at low tem- 
perature and decreases smoothly as the temperature is 
increased. This is also the case for the relaxation rate 
A3 of the fast-relaxing component. This smooth de- 
crease of these relaxation rates with temperature has 
been observed previously in magnetic materials^^'^^ 
and seems to roughly track the magnitude of the local 
field. It is possible that muon sites responsible for Ai 
and A3 lie further from the 2D planes than those sites 
giving rise to A2, and are thus less sensitive to 2D 
fluctuations, reducing the influence of any variation 
in r. The temperature evolution of Ai and A3 might 
then be expected to be dominated by the magnitude 
of A, which scales with the size of the local field and 
would therefore decrease as the magnetic transition is 
approached from below. 

The need for nonzero phases has been identified 
in previous studies of molecular magnets^'^'^^'^^, but 
never satisfactorily explained. One possible explana- 
tion for these might be that the muon experiences de- 
layed state formation. However, we can rule out the 
simplest model of this as the phases appear not to cor- 
relate with Vi. Such a correlation would be expected 
since a delay of to before entering the precessing state 
would give rise to a component of the relaxation func- 
tion ai{t) = cos [27rui {t + to)] = cos {27ruit + ^i), with 

(X which is not observed. This does not com- 
pletely rule out delayed state formation, as to could 
be a function of temperature (although this seems un- 
likely at these temperatures) . Nonzero phases are also 
sometimes observed when attempting to fit data with 
cosinusoidal relaxation functions from systems having 
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incommensurate magnetic structures. The phase then 
emerges as an artifact of fitting, as a cosine with a J 
phase shift approximates the zeroth-order Bessel func- 
tion of the first kind Jo{oJt) which is obtained from 
p+SR of an incommensurately-ordered system^^'^^. 
The Bessel function arises because the distribution of 
fields seen by muons at sites is asymmetric. However, 
attempts to fit the data with a pair of damped Bessel 
functions produced consistently worse fits than fits 
to Eq. (3), suggesting that a simple incommensurate 
structure is not a satisfactory explanation. It is also 
possible that several further magnetically-inequivalent 
muon sites exist, resulting in multiple, closely-spaced 
frequencies which give the spectra a more complex 
character which is not reflected in the fitting function. 
The simpler relaxation function would then obtain a 
better fit if the phase were allowed to vary. This has 
been observed^^, for example, in LiCr02. A final pos- 
sibility is that the distribution of fields at muon sites is 
asymmetric for another reason, perhaps arising from 
a complex magnetic structure. This may give rise to 
a Fourier transform which is only able to be fitted 
with phase-shifted cosines. However, the mechanism 
by which this would occur is unclear. 



B. Parametrizing exchange anisotropy 

The extent to which these systems approximate the 
2DSLQHA can be quantified by comparing the tran- 
sition temperature Tn to the exchange parameter J. 
The temperature Tn can be extracted using p+SR, 
whilst J can be obtained reliably from pulsed-field 
magnetization measurements^, heat capacity or mag- 
netic susceptibility. 

Mean-field theory predicts a simple relationship for 
the ratio of the transition temperature Tn and the 
exchange J given by^^ 



faTN 
J 



(5) 



where /cb is Boltzmann's constant, z is the number of 
nearest neighbors and S is the spin of the magnetic 
ions. In the pseudocubic [Cu(HF2)(pyz)2]X systems, 
S = ^ and z = 6, and Eq. (5) yields k^T^/J = 3. 
However, the reduced dimensionality increases the 
prevalence of quantum fluctuations, depressing the 
transition temperature and in [Cu(HF2)(pyz)2]BF4, 
we find k^T^/J ~ 0.25, which is indicative of large 
exchange anisotropy. 

Combining the experimental measures of Tn and J 
with the results of quantum Monte Carlo (QMC) sim- 
ulations allows us to deduce the exchange anisotropy 
Jx/J in the system Specifically, QMC simulations^'' 
for 2DSLQHA where 10"^ < J^/J <l are weh de- 
scribed by the expression 



J 



_ g6-47ip,/TN 



(6) 



Ps/T = 0.183 and h = 2.43. This expression allows a 
better estimate of k^T^/J in a 3D magnet: evaluat- 
ing for Jx/J = 1 yields /cbTn/J = 0.95. This is lower 
than the crude mean-field estimate because mean-field 
theory takes no account of fluctuations. Estimates of 
J for our materials, from pulsed magnetic field stud- 
ies except where noted, along with calculated Jx/J 
ratios, are shown in the summary tables throughout 
this paper. 

Another method of parametrizing the exchange 
anisotropy is to consider the predicted correlation 
length of two-dimensional correlations in the layers 
at the temperature at which we observe the onset 
of LRO. The larger this length, the better isolated 
the layers can be supposed to be. This can be es- 
timated by combining an analytic expression for the 
correlation length, in a pure 2DSLQHA^^, ^2D, with 
quantum Monte Carlo simulations to obtain an ex- 
pression^^'^^ appropriate for 1 < ^2d/<^ < 350,000, 



^ = 0.498ei-^^^^/^-^ 



1 - 0.44 



/cbT 
J 



O 



kBT 



where ps is the spin stiffness and 6 is a numerical con- 



stant. For S 



the appropriate parameters are 



J 

(7) 

where a is the square lattice constant, and T is the 
temperature. This formula yields ^2d(^n) ~ 0.5a for 
the mean-field model (/cbTn/J = 3), and ^2d(Tn) ^ a 
for /cbTn/J = 0.95 from quantum Monte Carlo simu- 
lations (i.e. Eq. (6) with J±/J = 1). By comparison, 
in [Cu(HF2)(pyz)2]BF4 Eq. (7) gives ^2d{Tn) « 50a, 
showing a dramatic increase in the size of correlated 
regions which build up in the quasi-2D layers before 
the onset of LRO. 



C. Nonmonotonic field dependence of Tn 

Although we expect the interplane exchange cou- 
pling J_i to have a large amount of control of the ther- 
modynamic properties of these materials, it may be 
the case that single-ion anisotropics are also respon- 
sible for deviations in the behaviour of our materials 
from the predictions of the 2DSLQHAF model. In 
particular, these anisotropics has been demonstrated 
to show a crossover to magnetic behaviour consistent 
with the 2D XY model^^. It was recently reported^^ 
that [Cu(HF2)(pyz)2]BF4 exhibits an unusual non- 
monotonic dependence of Tn as a function of applied 
magnetic field B [see Fig. 6(d)]. This behavior was ex- 
plained as resulting from the small XF-like anisotropy 
of the spin system in these systems. The physics of the 
unusual field-dependence then arises due to the dual 
effect of B on the spins, both suppressing the ampli- 
tude of the order parameter by polarizing the spins 
along a given direction, and also reducing the phase 
fluctuations by changing the order parameter phase 
space from a sphere to a circle. A more detailed ex- 
planation for the behavior^^ reveals that the energy 
scales of the physics are controlled by a Kosterlitz- 
Thouless-like mechanism, along with the interlayer ex- 
change interaction J±. 
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Figure 6. Sample TF p+SR data measured for [Cu(HF2)(pyz)2]BF4 in an applied field of 2 T are shown in (a) and 
(b). Data are shown in the 'rotating reference frame', rotating at 7]j x 1.9 T = 257 MHz, nearly cancelling out spin 
precession induced by the 2 T applied transverse field, (c) The evolution of the magnetic broadening a with T, showing 
a magnetic transition at 1.98 K in 2 T. (d) The B-T phase diagram from Ref. 12 showing the nonmonotonic behaviour 
at low applied magnetic field. In the key, HC is heat capacity, theory represents the results of computational modelling, 
and p^SR shows our results from TF measurements (see main text). 



The measurement of the B-T phase diagram in 
[Cu(HF2)(pyz)2]BF4 reported in Ref. 12 was made by 
observing a small anomaly in specific heat. In or- 
der to test whether the phase boundary could be de- 
termined using muons, we carried out transverse-field 
(TF) p+SR measurements using the LTF instrument 
at SpS. In these measurements, the field is applied per- 
pendicular to the initial muon spin direction, causing a 
precession of the muon-spins in the sum of the applied 
and internal field directed perpendicular to the muon- 
spin orientation. Example TF spectra measured in a 
field of 2 T are shown in Fig. 6 (a) and (b). We find 
that the spectra are well described by a function 

A{t) = A(0)e"^'^'/^ cos{27nyt + (^), (8) 

where the phase factor depends on the details of 
the detector geometry, and a is proportional to the 
second moment of the internal field distribution via 
cr^ = 7^(5^). Upon cooling through Tn we see a large 
increase in a, as shown in Fig. 6 (c). This approxi- 
mately resembles an order parameter, and we identify 
the discontinuity at the onset of the increase with Tn 
by fitting a with the above-T^ relaxation adding in 
quadrature to the additional relaxation present be- 
low the transition. The resulting point at T^{B = 
2 T) = 1.98(4) K is shown to be consistent with the 
predicted low- field phase boundary in Fig. 6 (d). A 
further point, identifiable by its vertical rather than 
horizontal error bar, was found by performing a field 
scan at a fixed temperature of T = 1.8 K. The field- 
dependence of the relaxation rate shows a sharp in- 
crease at the transition, at5 = 1.5±0.3T. 

Points derived from p+SR measurements possibly 
lie slightly lower in T than both that predicted by the- 
ory, and the line predicted on the basis of the specific 
heat measurements. The theoretical calculations use 



J/kB = 5.9 K and J_l/J = 2.5 x 10~^, whilst our esti- 
mates suggest J/kB = 6.3 K and Jj_/ J = 0.9 x 10~^. 
Performing these calculations for a purely 2D sys- 
tem results in the entire curve shifting to the left^^, 
and consequently the leftward shift of our data points 
is consistent with our finding of increased exchange 
anisotropy. It is clear that the TF p+SR technique 
may be used in future to measure the B-T phase dia- 
gram and enjoys some of the same advantages it has in 
ZF over specific heat and susceptibility in anisotropic 
systems. 



D. Muon response for T > Tn 

Above Tn, the character of the measured spectra 
changes considerably and we observe lower-frequency 
oscillations characteristic of the dipole-dipole interac- 
tion between muons and fluorine nuclei^^. The Cu^+ 
electronic moments, which dominate the spectra for 
T < Tn, are disordered in the paramagnetic regime 
and fluctuate very rapidly on the muon time scale. 
They are therefore motionally narrowed from the spec- 
tra, leaving the muon sensitive to the quasi-static nu- 
clear magnetic moments. 

A muon and nucleus interact via the two-spin 
Hamiltonian 

i>j 

where the spins Sij with gyromagnetic ratios jij are 
separated by the vector r. This gives rise to a preces- 
sion of the muon spin, and the muon-spin polarization 



9 



M 


X 


Tp-F (nm) 


PI (%) 


a (MHz) T (K) 


Cu 


BF4 


0.1038(1) 


77(1) 






Cu 


CIO4 


0.1081(2) 


70(1) 


yj.o t y± J 


0.1 


Cu 


PFe 


0.1044(2) 


74(2) 


0.29(2) 


5.2 


Cu 


AsFe 


0.1043(3) 


78(2) 


0.31(3) 


4.9 


Cu 


SbFe 


0.1047(2) 


64(1) 


0.30(1) 


5.0 


Cu 


BF4 


0.1042(1) 


76(1) 


0.26(1) 


26 


Cu 


CIO4 


0.1087(1) 


71(1) 


0.37(1) 


25 


Cu 


SbFe 


0.1080(3) 


59(1) 


0.26(1) 


30 


Cu 


NbFe 


0.1039(4) 


69(3) 


0.32(3) 


32 


Cu TaFe 


0.1039(2) 


78(2) 


0.26(3) 


32 


Ni 


SbFe 


0.1068(4) 


60(1) 


0.38(1) 


19 


Ni 


PFe 


0.1063(5) 


66(2) 


0.40(1) 


8.4 



Table II. Muon-fluorine dipole-dipole interaction fitted 
parameters in the family [M(HF2)(pyz)2]X, extracted 
from fitting data to Eq. (12). In addition to separation 
by metal ion, Cu compounds are grouped by the temper- 
ature at which the measurement was made: those com- 
pounds measured over a range of temperatures appear in 
both sections of the table. 




5 10 



t (ps) 

Figure 7. Data taken at T = 15 K > Tn = 1.4 K for 
[Cu(HF2)(pyz)2]BF4, showing Fp oscillations along with 
a fit to Eq. (11). The inset shows the energy levels present 
in a simple system of two 5* = | spins, along with the 
allowed transitions. 



along a quantization axis z varies with time as 



(10) 



where N is the number of spin states, |m) and \n) are 
eigenstates of the total Hamiltonian cFq is the Pauli 
spin matrix corresponding to the direction and ()^ 
represents an appropriately- weighted powder average. 
The vibrational frequency of the muon-fluorine bond 
exceeds by orders of magnitude both the frequencies 
observable in a p+SR experiment, and the frequency 
appropriate to the dipolar coupling in Eq. (9); the 
bond length probed via these entangled states is thus 
time- averaged over thermal fluctuations. Fluorine is 
an especially strong candidate for this type of interac- 
tion firstly because it is highly electronegative causing 
the positive muon to stop close to fluorine ions, and 
secondly because its nuclei are 100% ^^F, which has 

i = \. 

Data were fitted to a relaxation function 

A{t) = Ao(pie-^---*D,(t) +p2e-"'*') + Age'^^^*, 

(11) 

where the amplitude fraction pi ^ 70% reflects the 
muons stopping in a site or set of sites near to a 
fluorine nucleus, which result in the observed oscil- 
lations Dz(t)] the weak relaxation of the muon spins 
is crudely modelled by a decaying exponential. The 
fraction p2 ^ 30% describes those muons stopping in 
a class of sites primarily influenced by the randomly- 
orientated fields from other nuclear moments, giving 
rise to a Gaussian relaxation with a ^ 0.'^ MHz. Ex- 
ample data and a fit are shown in Fig. 7, whilst param- 
eters extracted by fitting this function to data from 
each compound are shown in Table II. 

Fits to a variety of different Dz{t) functions were 
attempted, including that resulting from a simple Fp 



bond (previously observed in some polymers^^) and 
the better-known FpF complex comprising a muon 
and two fluorine nuclei in linear symmetric configu- 
ration, which is seen in many alkali fluorides^^. This 
latter model was also modified to include the possi- 
bilities of asymmetric and nonlinear bonds. Previ- 
ous measurements^^ made in the paramagnetic regime 
of [Cu(HF2)(pyz)2]C104 suggested that the muon 
stopped close to a single fluorine in the HF2 group and 
also interacted with the more distant proton. This in- 
teraction is dominated by the F-p coupling and, for 
our fitting, the observed muon-fluorine dipole-dipole 
oscillations were found to be well described by a sin- 
gle Fp interaction damped by a phenomenological re- 
laxation factor. For such Fp entanglement, the time 
evolution of the polarization is described by 



(12) 



where ixi = 2, 2x2 = 1 and = 2. The frequen- 
cies ujj = jujdf^^ where uj^ = /io7ti7F^/47rr^, in which 
7f = 27rx 2.518x10^ MHzT"^ is the gyromagnetic ra- 
tio of a ^^F nucleus^^, and r is the muon-fluorine sep- 
aration. These three frequencies arise from the three 
transitions between the three energy levels present in 
a system of two entangled S = ^ particles (see in- 
set to Fig. 7). The fact that the relaxation func- 
tion is similar in all materials in the series, includ- 
ing [Cu(HF2)(pyz)2]C104 which is the only compound 
studied without fluorine in its anion, (the only differ- 
ence being a slight lengthening of the p-F bond, and 
with no significant change in oscillating fraction) sug- 
gests that the muon site giving rise to the Fp oscilla- 
tions in all systems is near the HF2 bridging ligand. 

The temperature evolution of the Fp signal was 
studied for T < 300 K in [Cu(HF2)(pyz)2]BF4 and 
[Cu(HF2)(pyz)2]C104. In both cases, the dipole- 
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Figure 8. (a) Muon-fluorine dipole-dipole oscillations in [Cu(HF2)(pyz)2]BF4 over a temperature range 2 < T < 300 K. 
Asymmetry spectra are displaced vertically so as to approximately align with the temperature scale on plots (b) and (c). 
Ticks on the ^/-axis of (a) denote 1% asymmetry. Plot (b) shows the fitted value for rp-p as a function of temperature. 
The line shown is a fit to the low-T points with a scaling law, Eq. (14). The upper x-axis shows values of the 
dipole frequency, z/d, which correspond to the lower x-axis values of rp-p. Plot (c) shows fitted relaxation rates Ap-p 
and Ao, referring to the relaxation of the Fp function Dzit) in Eq. (11), and the pure relaxation in Eq. (13). Shaded 
regions indicate temperatures where Q = 2n\Y-^/uJd > 2, roughly parametrizing the disappearance of the oscillations, 
(d), (e) and (f) follow (a), (b) and (c), but show data for [Cu(HF2)(pyz)2]C104. (g), (h) and (i) similarly, but for 
[Cu(HF2)(pyz)2]SbF6 over the range 5 < T < 250 K. The 5 K A(t) plot is omitted because the background is raised 
substantially by approach to the transition to LRO. 
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material ro (nm) a (10 ^nmK ^) 

[Cu(HF2)(pyz)2]BF4 0.10376(3) 3.96(6) 
[Cu(HF2)(pyz)2]C104 0.10842 2.7212 

PVDF^2 0.10914 1.9488 



Table III. Fitted values obtained by fitting muon-fluorine 
bond lengths with a scaling law, as Eq. (14). 

dipole oscillations disappear gradually in a tempera- 
ture range 150 < T < 250 K, with oscillations totally 
absent in the center of this range, followed by reap- 
pearing as temperature is increased further. Plots of 
A{t) spectra at a variety of temperatures are shown 
in Fig. 8 (a) and (d). The data were initially fit- 
ted to Eq. (11), with all parameters left free to vary. 
The temperature-evolution of the muon-fluorine bond 
length, rp can be seen in Fig. 8 (b) and (e). The 
spectra were also fitted with 

A{t) = Ao(pie-^°^ +p2e-"''') + Abge-^^^^ (13) 

a sum of an exponential and a Gaussian relaxation, 
which might be expected to describe the data in the 
region where the oscillations vanish. Both this relax- 
ation and that extracted from Eq. (11) are plotted in 
Fig. 8 (c) and (f), labelled Aq and Ap p respectively. 

This bond length appears to grow and then shrink 
by nearly 20% over the 100 K range where the os- 
cillations fade from the spectra and reappear. This 
variation is significantly larger than any variation in 
crystal lattice parameters which would be expected. 
Since the oscillations visibly disappear from the mea- 
sured spectra, results from fitting with an oscillatory 
relaxation function are artifacts of the fitting proce- 
dure: since the frequencies scale with 1/r^, increasing 
bond length together with the associated relaxation 
rate fits the data with a suppressed oscillatory signal. 
This can be approximately quantified by examining 
the ratio Q = 27rAF p/c^d, where a large value indi- 
cates that the function relaxes significantly before a 
single Fp oscillation is completed. The shaded re- 
gions in Fig. 8 show where Q > 2, which acts as an 
approximate bound on where the parametrization in 
Eq. (11) would be expected to fail. In the low-T re- 
gion where Q < 2, the bond lengths appear to scale 
roughly as T^, which has previously been observed in 
fluoropolymers^^ . Parameters extracted from fitting 
to 

rp ^ = aT2+ro (14) 

are shown in Table III. 

The observation in these two samples of Fp oscilla- 
tions which disappear and then reappear is puzzling. 
While we have not identified a definitive mechanism, 
we can probably rule out an electronically mediated 
effect since, for T ^ T^^ the Cu moment fluctuations 
will be outside the muon time-window. An explana- 
tion could involve nearby nuclear moments, possibly 
influenced by a thermally-driven structural distortion 
or instability. 



A similar study of [Cu(HF2)(pyz)2]SbF6 is shown 
in Fig. 8 (g), (h) and (i). In this material, the os- 
cillations appear not to vanish over the temperature 
range studied, though we cannot rule out a brief dis- 
appearance Sit T ^ 200 K. Instead, the oscillations 
show an apparently monotonic increase in damping 
with temperature, and the fitted bond length does 
not follow Eq. (14). The shaded region in Fig. 8 (h) 
and (i) has no upper bound, though we cannot rule 
out a constraint at T > 250 K. The pure relaxation 
Ao is omitted because there is no region where the Fp 
oscillations are sufficiently damped for Eq. (13) to be 
a good parametrization. 



E. Muon site determination 

Combining the data measured above and below the 
transition in these materials allows us to attempt to 
construct a self-consistent picture of possible muon 
sites. The observed dipole-dipole observations above 
Tn suggest that at least one muon stopping site is 
near a fluorine ion. We consider three classes of prob- 
able muon site: Class I sites near the fluorine ions in 
the HF2 groups. Class II sites near the pyrazine rings, 
and Class III sites near the anions at the centre of 
the pseudocubic pores. Comparison of Tables I and 
II show that the dominant amplitude component for 
T > Tn arises from dipole-dipole oscillatory compo- 
nent pie~^^Dz{t) and from the fast-relaxing compo- 
nent pse~^^^ for T < Tn, and that these are compa- 
rable in amplitude. It is plausible therefore, to sug- 
gest that these two signals correspond to contribu- 
tions from the same Class I muon sites near the HF2 
groups. Moreover, the analysis of the T > Tn spec- 
tra in the previous section implies that this site lies 

F ^0.11 nm from an F in the HF2 groups. The 
remainder of the signal (the oscillating fraction below 
Tn and the Gaussian relaxation above) can also be 
identified, suggesting that the sites uncoupled from 
fluorine nuclei (Classes II and/or III) result in the 
magnetic oscillations observed for T < Tn- 

We note further that the evidence from Fp oscil- 
lations makes the occurence of Class III muon sites 
unlikely. The fact that spectra observed for T > Tn 
in the X~ = CIO4 material are nearly identical to 
those in all other compounds, in which X contains flu- 
orine, suggest that the muons do not stop near the an- 
ions. Moreover, as discussed in Sec. IV and V below, 
we observe no Fp oscillations in [Cu(pyz)2(pyo)2]X2 
{Y- = BF4, PFg) or [Cu(pyo)6](BF4)2, suggest- 
ing that muons do not stop preferentially near these 
fluorine-rich anions either. We therefore rule out the 
existence of Class III muon sites and propose that the 
magnetic oscillations measured for T < Tn most prob- 
ably arise due to Class II sites found near the pyrazine 
ligands. 

Below Tn, the measured muon precession frequen- 
cies allow us to determine the magnetic field at these 
Class II muon sites via u = 7p5/27r. Simulating the 
magnetic field inside the crystal therefore allows us to 
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Figure 9. Probability density functions of muon precession 
frequencies at positions close to likely muon stopping sites 
in [Cu(HF2)(pyz)2]X with Cu^^ moments /x = /xb- The 
graphs show dipole fields near the fluorine or oxygen atoms 
in the negative anions; near the fluorine atoms in the bi- 
fluoride ligands; and near the carbon and nitrogen atoms, 
as a proxy for proximity to the pyrazine ring. The type of 
line indicates the compound for which the calculation was 
performed. The muon site is constrained to be close to 
particular atoms, indicated by line color. The shaded ar- 
eas indicate ranges of fitted frequencies as T ^ 0; two fre- 
quencies z^^2 ^ represent those observed where X~ — BF4 , 
and those observed in the X~ — PF^ analogue. 



compare these 5-fields with those predicted for likely 
magnetic structures and may permit us to constrain 
the ordered moment. For the case of our ZF measure- 
ments in the antiferromagnetic state, the local mag- 
netic field at the muon site .Biocai is given by 



local 



dipole 



B 



hyperfine •> 



(15) 



where -Bdipoie is the dipolar field from magnetic ions 
located within a large sphere centred on the muon 
site and ^hyperfine the contact hyperfine field caused 
by any spin density overlapping with the muon wave- 
function. This spin density is difficult to estimate ac- 
curately, particularly in complex molecular systems, 
but it is probable for insulating materials such as these 
that the spin density on the copper ion is well localised 
and so we ignore the hyperfine contribution in our 
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Figure 10. Probability density functions for muons in 
the putative oscillating sites near the pyrazine rings in 
[Cu(HF2)(pyz)2]BF4 (a) for muon precession frequency v 
assuming that the moment on the copper site \i — 
created fom a histogram of dipole fields evaluated using 
Eq. (16) at points satisfying the constraints detailed in 
the text; and (b) for moment on the copper sites given 
the frequencies actually observed, evaluated using the pdfs 
in (a) and Eq. (18). Lines represent trial magnetic struc- 
tures. All exhibit antiferromagnetic coupling through both 
the bifluoride and pyrazine exchange paths, whilst (i) has 
copper moments pointing at 45° to the pyrazine grid, (ii) 
has moments along one of the pyrazine grid directions (a 
or 6), and (iii) has copper moments pointing along the 
bifluoride axis (c). 



analysis. The dipole field ^dipole is a function of the 
coordinate of the muon site r^, and comprises a vec- 
tor sum of the fields from each of the magnetic ions 
in the crystal approximated as a point dipole, so that 



^dipoie(r'n) = M [3(Ai • r)r 



(16) 



where r = — is the relative position of the muon 
and the i^^ ion with magnetic moment fii = iifii^ and 
i is an index implying summation over all of the ions 
which make up the crystal. 

Although these materials are known to be anti- 
ferromagnetic from their negative Curie-Weiss tem- 
peratures and zero spontaneous magnetization at low 
temperatures^^'^^, their magnetic structures are un- 
known. Dipole field simulations were therefore per- 
formed for a variety of trial magnetic structures with 
jii = jiiB- We analyse the results of these calculations 
using a probabilistic method. We begin by allow- 
ing the possibility that the magnetic precession sig- 
nal could arise from any of the possible classes of 
muon site identified above. Random positions in the 
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unit cell were generated and dipole fields calculated 
at these. To prevent candidate sites lying too close to 
atoms we constrain all sites such that r^, a > 0.1 nm 
where A is any atom. Possible Class I muon sites 
were identified with p = ±0.01 nm (where tq is 
the muon-fluorine distance established from Fp os- 
cillations) and possible Class II sites were selected 
with the constraint that 0.10 < c,n ^ 0.12 nm. 
The predicted probability density function (pdf) of 
muon precession frequencies (resulting from the mag- 
nitudes of the calculated fields) are plotted in Fig. 9, 
with the observed frequencies superimposed. Results 
are shown for a trial magnetic structure compris- 
ing copper spins lying in the plane of the pyrazine 
layers and at 45° to the directions of the pyrazine 
chains, and with spins arranged antiferromagnetically 
both along those chains and along the HF2 groups. 
This candidate structure is motivated by analogy with 
[Cu(pyz)2](C104)2, which also comprises Cu^+ ions in 
layers of 2D pyrazine lattices^^, and with the parent 
phases of the cuprate superconductors^^, which are 
also two-dimensional Heisenberg systems of S = ^ 
Cu^+ ions. Other magnetic structures investigated 
give qualitatively similar results. From Fig. 9 it is 
clear that the only sites with significant probability 
density near to the observed frequencies are those ly- 
ing near the anions (i.e. Class III sites) which we have 
argued are not compatible with our data. The more 
plausible muon sites correspond to higher frequencies 
than those observed. Our conclusion is that it is likely 
that the Cu^+ moments are rather smaller than the 
/iB assumed in the initial calculation. 

If we accept that the muon sites giving rise to mag- 
netic precession are near the pyrazine groups then we 
may use this calculation to constrain the size of the 
copper moment. Since u is obtained from experiment, 
what we would like to know is ^(/i|i/), the pdf of copper 
moment /i given the observed u. This can be obtained 
from our calculated /{u/ja) using Bayes' theorem^^, 
which yields 



(17) 



where we have assumed a prior probability for the cop- 
per moment that is uniform between zero and /imax- 
We take /imax = 2/iB, although our results are in- 
sensitive to the precise value of /imax as long as it 
is reasonably large. When multiple frequencies {ui} 
are present in the spectra, it is necessary to multiply 
their probabilities of observation in order to obtain 
the chance of their simultaneous observation, so we 
evaluate 



g{fi\{iyi}) (X 



n/ 



f{Pi/^)dPi, (18) 



where Aui is the error on the fitted frequency. Re- 
sults are shown in Fig. 10, along with the dipole field 
pdfs which gave rise to them. By inspection of the 
pdfs, the copper moment is likely to be /i < 0.5/iB- 




Figure 11. Structure of [Cu(pyz)2(pyo)2](BF4)2- Cop- 
per ions lie in 2D square layers, bound by pyrazine rings. 
Pyridine- A/'-oxide ligands protrude from the coppers in 
a direction approximately perpendicular to these layers. 
Tetrafluoroboride ions fill the pores remaining in the struc- 
ture. Ion sizes are schematic; copper ions are shown twice 
as large for emphasis, and hydrogens have been omitted 
for clarity. 



The dipole field simulations results also lend weight 
to our contention that the oscillatory signal cannot 
arise from the sites that also lead to the Fp compo- 
nent above-TN- If this were the case then the most 
likely moment on the copper would be /icu ^ 0.2/iB, 
which seems unreasonably small. We note that mo- 
ment sizes of /i < 0.5/iB were also observed for the 
2DSLQHA system La2Cu04 (a recent estimate^^ from 
neutron diffraction gave [0.42 ± 0.01] /ib), despite the 
predictions of 0.6/iB from spin wave theory and Quan- 
tum Monte Carlo^^. It was suggested in that case^^ 
that disorder might play a role in reducing the mo- 
ment sizes; an additional possible mechanism for this 
suppression is ring exchange^^'^^. 

One limitation of this analysis is that the mecha- 
nism for magnetic coupling of copper ions through the 
pyrazine rings is postulated to be via spin exchange, 
in which small magnetic polarisations are induced on 
intervening atoms^^. Density functional theory calcu- 
lations estimate that these are small, with the nitro- 
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Figure 12. Example data and fits for [Cu(pyz)2(pyo)2] 1^2- From left to right: (a) and (d) show sample asymmetry 
spectra A{t) for T < Tn along with a fit to Eq. (19); (b) and (e) show the frequency zy as a function of temperature; and 
(c) and (f) show relaxation rates as a function of temperature. The relaxation rates Ai, associated with the oscillation, 
and As, the fast-relaxing initial component, do not vary significantly with T and are not shown. Only a slight trend in 
A2 in the Y~ = PF^ material [graph (f)] is observed. 



gen and carbon moments estimated at /ic ^ 0.01/iB 
and /iN ^ 0.07/iB, respectively^^. However, their ef- 
fect may be non-negligible: they may be significantly 
closer to the muon site than a copper moment, and 
dipole fields fall off rapidly, as 1/r^. Further, since 
much of the electron density in a pyrazine ring is de- 
localised in Tt-orbitals, the moments may not be point- 
like, as assumed in our dipole field calculations. Fur- 
ther, this may lead to overlap of spin density at the 
muon site and result in a nonzero hyperfine field. 



IV. [Cu(pyz)2(pyo)2]y2 

In this section, we report the magnetic behavior 
of another family of molecular systems which shows 
quasi-2D magnetism, but for which the interlayer 
groups are very different and arranged in a com- 
pletely different structure, resulting in a 2D coordi- 
nation polymer. This system is [Cu(pyz)2(pyo)2] 1^2, 
where Y~ = BF4 , PFg . As with the previous case, 
S = ^ Cu^+ ions are bound in a 2D square lattice of 
[Cu(pyz)2]^^ sheets lying in the a6-plane. Pyridine- TV- 
oxide (pyo) ligands [shown in Fig. 1 (b)] protrude from 
the copper ions along the c-direction, perpendicular 
to the a6-plane in the Y~ = PF^ material, but mak- 
ing an angle /3 — 90 ^ 29° with the normal in Y~ = 
BF4 . The anions then fill the pores remaining in the 
structure. The structure of [Cu(pyz)2(pyo)2](BF4)2 is 
shown in Fig. 11. 

In a typical synthesis, an aqueous solution of Cu Y2 
hydrate ( Y~ = BF J or PFg ) was combined with an 
ethanol solution that contained a mixture of pyrazine 



and pyridine- A/"- oxide or 4-phenylpyridine-A/'-oxide. 
Deep blue-green solutions were obtained in each case, 
and when allowed to slowly evaporate at room temper- 
ature for a few weeks, dark green plates were recovered 
in high yield. Crystal quality could be improved by se- 
quential dilution and collection of multiple batches of 
crystals from the original mother liquor. The relative 
amounts of pyz and pyo were optimized in order to 
prevent formation of compounds such as Cuy2(pyz)2 
or [Cu(pyo)6] Y2. 

Samples were measured in the LTF apparatus at 
SpS. Example data measured on [Cu(pyz)2(pyo)2] Y2 
are shown in Fig. 12, where we observe oscillations in 
A{t) at a single frequency below T^. Data were fitted 
to a relaxation function 



A{t) = Ao {pi cos{2nuit)e-^'* ^p2e-^^* pse 



(19) 



The small amplitude fraction pi < 10% for both sam- 
ples refers to muons stopping in a site or set of sites 
with a narrow distribution of quasi-static local mag- 
netic fields, giving rise to the oscillations; p2 ~ 50% is 
the fraction of muons stopping in a class of sites giving 
rise to a large relaxation rate 30 < A < 60 MHz and 
Ps ~ 50% represents the fraction of muons stopping 
in sites with a small relaxation rate A3 ~ 1 MHz. The 
data from these compounds fit best with ^ = 0, and 
it is thus omitted from this expression. Frequencies 
obtained from fitting the data to Eq. (19) were then 
modelled with Eq. (4). The results of these fits are 
shown Fig. 12, and Table IV. 

Our results show that [Cu(pyz)2(pyo)2](BF4)2 has 
a transition temperature Tn = 1.5dz0.1K and a 
quasi-static magnetic field at the muon site ui{T = 



0) = 1.4 ± 0.1 MHz. No quantities other than ui 
show a significant trend in the temperature region 
0.1 < T < 1.6 K. Above the transition, purely re- 
laxing spectra are observed, displaying no Fp oscilla- 
tions. As suggested above, this makes the existence of 
muon sites near the anions unlikely. We find a critical 
exponent of (3 = 0.25 ± 0.10, where the large uncer- 
tainty results in part from the difficulty in fitting the 
A{t) data in the critical region. 

Our results for [Cu(pyz)2(pyo)2](PF6)2 show that 
the transition temperature is slightly higher at Tn = 
1.72 ± 0.02 K and the oscillations occur at a lower 
frequency of ui{T = 0) = 1.07 ± 0.03 MHz. The re- 
laxation rates A2 and A3 also decrease in magnitude as 
temperature is increased, settling on roughly constant 
values A2 ^ 0.6 MHz and A3 ^ 15 MHz for T > Tn- 
No other quantities show a significant trend in the 
temperature region 0.2 < T < 1.7K. Above the tran- 
sition, relaxing spectra devoid of Fp oscillations are 
again observed. The critical exponent (3 = 0.22 ±0.02. 

The small amplitude of the oscillations, common 
to both samples, might be explained in a number of 
ways. The materials may undergo long-range ordering 
but there may be an increased likelihood of stopping 
in sites where the magnetic field nearly precisely can- 
cels. Alternatively, a range of similar muon sites may 
be present with a large distribution of frequencies, or 
alternatively the presence of dynamics, washing out 
any clear oscillations in large fractions of the spectra 
and instead resulting in a relaxation. Finally, we can- 
not exclude the possibility that only a small volume 
of the sample undergoes a magnetic transition; this 
may indicate the presence of a small impurity phase, 
possibly located at either grain boundaries or, given 
that this is a powder sample, near the crystallites' 
surfaces. We note also that the behavior of fitted pa- 
rameters in these materials is qualitatively similar to 
that reported in CuCl2(pyz), where there is also a rel- 
atively small precessing fraction of muons and little 
variation in relaxation rates as Tn is approached from 
below^^. 
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Figure 13. Structure of [Cu(pyo)6](BF4)2, viewed along 
the three-fold (c-) axis, after Ref. 44. Copper ions are 
surrounded by octahedra of six oxygens, each part of a 
pyridine- iV-oxide ligand; [Cu(pyo)6]^^ complexes space- 
pack with BF^ stabilising the structure. Ion sizes are 
schematic; copper ions are shown twice as large for em- 
phasis, and hydrogens have been omitted for clarity. As 
indicated in the top left, the a and b directions lie in the 
plane of the paper, separated by 7 = 120°, whilst the c 
direction is out of the page. 



V. [Cu(pyo)3]Z2 

The next example is not a coordination polymer, 
but instead forms a three-dimensional structure of 
packed molecular groups. The molecular magnet 
[Cu(pyo)6]^2, where Z' = BF4 , CIO3", PFg , com- 
prises Cu^+ ions on a slightly distorted cubic lattice, 
located in [Cu(pyo)6]^^ complexes, and surrounded by 
octahedra of oxygen atoms^^. The structure is shown 
in 13. This approximately cubic structure, which 
arises from the molecules' packing, might suggest that 
a three-dimensional model of magnetism would be ap- 
propriate. In fact, although the observed bulk prop- 
erties of [M(pyo)6]X2 where M'^^ = Co^+, Ni^^ or 
Fe^+ are largely isotropic, but the copper analogues 
display quasi-low-dimensional, S = ^ Heisenberg an- 
tiferromagnetism^^. Weakening of superexchange in 
certain directions, and thus the lowering of the sys- 



tems' effective dimensionality, is attributed to length- 
ening of the superexchange pathways resulting from 
Jahn-Teller distortion of the Cu-0 octahedra, which 
is observed in structural and EPR measurements^^'^^. 
At high temperatures, (T > 100 K), these distortions 
are expected to be dynamic but, as T is reduced (to 
^ 50 K), they freeze out. The anion Z~ determines 
the nature of the static Jahn-Teller elongation. The 
Z~ = BF4 material displays ferrodistortive order- 
ing which, in combination with the antiferromagnetic 
exchange, gives rise to 2D Heisenberg antiferromag- 
netic behavior"^^'"^^. By contrast, Z' = CIO4 , NO3 
(neither of which is investigated here) display antifer- 
rodistortive ordering^^, which gives rise to quasi- ID 
Heisenberg antiferromagnetism^^ . All of the samples 
investigated were measured in the LTF spectrometer 
at SpS. 
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Figure 14. Example data and fits for [Cu(pyo)6](BF4)2. From left to right: (a) shows sample asymmetry spectra A{t) 
for T < Tn and T > Tn, along with fits to Eq. (19) and Eq. (20), respectively; (b) shows frequency as a function 
of temperature; and (c) shows relaxation rate A2 as a function of temperature; Ai was held fixed during the fitting 
procedure, and A3 persists for T > Tn- In the z^(T) plot, error bars are included on the points but in most cases they 
are smaller than the marker being used. 



In the Z~ = BF4 compound, below a temperature 
Tn, a single oscillating frequency is observed, indi- 
cating a transition to a state of long-range magnetic 
order. Example data above and below the transition, 
along with fits, are shown in Fig. 14 (a). Data were 
fitted to Eq. (19), and the frequencies extracted from 
the procedure fitted as a function of temperature to 
Eq. (4)as shown in Fig. 14 (b). This procedure iden- 
tifies a transition temperature Tn = 0.649 ± 0.005 K. 
The only other parameter found to vary significantly 
in the range 20 mK < T < Tn was A2, shown 
in as shown in Fig. 14 (c). The fitted parameters 
are shown in Table IV. We may compare this with 
the result of an earlier low-temperature specific heat 
study^^ which found a very small A-point anomaly at 
Tn = 0.62±0.01 K, slightly lower than our result. Fit- 
ting the magnetic component of the heat capacity with 
the predictions from a two-dimensional Heisenberg 
antiferromagnet gives J//cb = — 1.10 ± 0.02 K, and 
similar analysis of the magnetic susceptibility^^ yields 
J/kB = — 1.08 ± 0.03 K. Using these values, together 
with the muon estimate of Tn and Eq. (6), allows us to 
estimate the inter-plane coupling, J±/ J = 0.26±0.01. 
(Using the value of Tn from heat capacity results in 
an estimate J±/J = 0.21 ± 0.02.) 

For temperatures Tn < T < 1 K, the spectra are 
well described by a relaxation function 

A{t) = Ao(pie-^^ +p2e-"''') + Abg, (20) 

comprising an initial fast-relaxing component with 
A ~ 20 MHz, and a Gaussian relaxation with a ^ 
0.4 MHz corresponding to the slow depolarisation of 
muon spins due to randomly-orientated nuclear mo- 
ments. 

The Z~ = ClO^ material also shows evidence for a 
magnetic transition, although in this case we do not 
observe oscillations in the muon asymmetry. Instead 
we measure a discontinuous change in the relaxation 
which seems to point towards an ordering transition. 
Example asymmetry spectra are shown in Fig. 15 (a) 
and data at all measured temperatures are well de- 



scribed with the relaxation function 

A{t) = Aoe-^' + Abg . (21) 

Evidence for a magnetic transition comes from the 
temperature evolution of A [Fig. 15(b)], where we 
see that the relaxation decreases with increasing tem- 
perature until it settles at T ^ 0.3 K, on a value 
A 0.5 MHz. It is likely that this tracks the internal 
magnetic field inside the material, and is suggestive of 
Tn = 0.30 ±0.01 K. 

The final member of this family studied, Z~ = PF^, 
shows no evidence for a magnetic transition over the 
range of temperatures studied, 0.02 K < T < 1 K. An 
example spectrum is shown in Fig. 16. The data re- 
semble the above-transition data measured in the BF4 
and CIO3 compounds and it therefore seems likely 
that the paramagnetic state persists to the lowest tem- 
perature measured. 



VI. Ag(pyz)2(S208) 

The examples so far have used Cu^+ (3d^) as the 
magnetic species. An alternative strategy is to employ 
Ag^+ (4d^) which also carries an /S = ^ moment. This 
idea has led to the synthesis of Ag(pyz)2(S2 08), which 
comprises square sheets of [Ag(pyz)2]^^ units spaced 
with 8208^" anions"^^. Each silver ion lies at the cen- 
tre of an elongated (AgN402) octahedron, where the 
Ag-N bonds are significantly shorter than the Ag-0. 
Preparation details can be found in Ref. 49. 

Measurements were made using the GPS instru- 
ment at SpS. Example muon data, along with fits to 
various parameters, are shown in Fig. 17. Asymmetry 
oscillations are visible in spectra taken below a tran- 
sition temperature Tn- The data were fitted with a 
relaxation function 

A{t) = Ao {pi cos{2nuit + ^i)e"^^^ + P2e~^^^) + Abg, 

(22) 
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Figure 15. Example data and fits for [Cu(pyo)6] (€103)2 • Representative A{t) spectra for T < Tn and T > Tn, along 
with fits to Eq. (21), are shown in (a), whilst the value of the relaxation rate, A, as a function of temperature is shown 
in (b). 



comprising a single damped oscillatory component, a 
slow-relaxing component, and a static background sig- 
nal. The onset of increased relaxation Ai as the tran- 
sition is approached from below leads to large statisti- 
cal errors on fitted values, as is evident in Fig. 17 (c). 
The relaxation A2 decreases with increasing tempera- 
ture. Fitting to Eq. (4) allows the critical parameters 
P = 0.19±0.02 and Tn = 7.8±0.3 K to be determined. 
Fitted values are shown in Table IV. 

The in-plane exchange J is too large to be deter- 
mined with pulsed fields^^: M{B) does not saturate 
in fields up to 64 T. However, fitting x{T) data al- 
lows an estimate of the exchange J//cb 53 K (and 
thus, in conjunction with g measured by EPR, the 
saturation field Be is estimated to be 160 T). Thus, 
kBT^/J = 0.148 ± 0.006. Estimation of the exhange 
anisotropy with Eq. (6) yields \ J_l/J\ ~ 10~^, but this 
very small ratio of ordering temperature to exchange 
strength is outside the range in which the equation 
is known to yield accurate results. The alternative 
method of parametrizing the low dimensionality in 
terms of correlation length at the Neel temperature 
(see Sec. IIIB) yields ^{T^)/a = 1000 ± 300. 
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Figure 16. Sample asymmetry spectrum for 

[Cu(pyo)6](PF6)2 measured at T = 0.02 K. The 
spectra remain indistinguishable from this across the 
range of temperatures examined, 0.02 K < T < 1 K. 



VII. [Ni(HF2)(pyz)2]X 

In order to investigate the influence of a dif- 
ferent spin state on the magnetic cation in the 
[M(HF2)(pyz)2]X architecture the [Ni(HF2)(pyz)2]X 
{X~ = PFg , SbFg ) system has been synthesised^^. 
These materials are isostructural with the copper fam- 
ily discussed in Sec. Ill, but contain 5* = 1 Ni^+ 
cations. 

Data were taken using the GPS spectrometer at 
PSI. Example data are shown in Fig. 18. We observe 
oscillations at two frequencies below the materials' re- 
spective ordering temperatures. Data were fitted with 
a relaxation function 

A{t) = Ao [pie~^'* cos{27TUit) + j92e"^'^ cos(27rP2i^i^) 

+P3e"^^'] + Abge-^b^^ (23) 

Of those muons which stop in the sample, pi ^ 
25% indicates the fraction of the signal correspond- 
ing to the low- frequency oscillating state with ui (T = 
0) ^ 12.3 MHz; p2 ^ 10% corresponds to muons 
stopping in the high-frequency oscillating state with 
V2{T = 0) ^ 9.0 MHz; and ps ~ 65% represents 
muons stopping in a site with a large relaxation rate 
XsiT = 0) ^ 70 MHz. The frequencies were ob- 
served to scale with one-another, and consequently 
the second frequency was held in fixed proportion 
1^2 = P2^i during the fitting procedure. The only 
other parameter which changes significantly in value 
below Tn is A3, which decreases with a trend qualita- 
tively similar to that of the frequencies. Fitting the 
extracted frequencies to Eq. (4) allows the transition 
temperature Tn = 12. 25 ±0.03 K and critical exponent 
(3 = 0.34 ± 0.04 to be extracted. In contrast to the 
copper family studied in Sec. HI, the relation X (x 
holds true, suggesting that a field distribution whose 
width diminishes with increasing temperature is re- 
sponsible for the variation in A, and that dynamics 
are relatively unimportant in determining the muon 
response. This is shown graphically in the inset to 
Fig. 18 (f), where a plot of frequency against relax- 
ation rate lies on top of a line representing a A = 
relationship. The phase <j) required in previous fits 
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Figure 17. Example data and fits for Ag(pyz)2(S208). From left to right: (a) shows sample asymmetry spectra A{t) 
for T < Tn along with a fit to Eq. (22); (b) shows the frequency as a function of temperature; and (c) shows relaxation 
rates as a function of temperature. Filled circles show the relaxation rate Ai, which relaxes the oscillation. The filled 
triangles correspond to the relaxation A2. 



material z/i (MHz) pi p2 P3 Tn (K) /3 a J/fe (K) \J±/J\ 

[Cu(pyz)2(pyo)2](BF4)2 1.4(1) < 10 50 50 1.5(1) 0.25(10) 1.6(3) 

[Cu(pyz)2(pyo)2](PF6)2 0.64(1) < 10 50 50 1.72(2) 0.22(2) 1.1(3) 8.1(3)* 2 x 10"^ 

[Cu(pyo)6](BF4)2 0.599(5) 15 30 55 0.649(5) 0.29(1) 1.7(1) 1.09(2) ^ 0.23(2) 

Ag(pyz)2(S208) 2.45(4) 30 70 - 7.8(3) 0.19(2) 3(2) 52.7(3) - 10"^ 



Table IV. Fitted parameters for molecular magnets in the family [Cu(pyz)2(pyo)2] 1^2, [Cu(pyo)6](BF4)2 and 
Ag(pyz) 2(8208). The first parameters shown relate to fits to Eq. (19) (for the first three rows) or Eq. (22) [for 
Ag(pyz)2(S208)]. This allows us to derive frequencies at T = 0, Ui; and probabilities of stopping in the various classes 
of stopping site, p^, in percent. Then, the temperature dependence of z/^ is fitted with Eq. (4), extracting values for 
the Neel temperature, Tn, critical exponent /3 and parameter a. Finally, the quoted J//cb is obtained from pulsed- field 
experiments^ The asterisk (*) indicates the value of J was obtained using gab = 2.04(1) in conjunction with single-crystal 
pulsed-field data. The dagger (^) indicates that the value of J is extracted from heat capacity and susceptibility from 
Ref. 45.). The ratio of inter- to in-plane coupling, J±/ J, is obtained by combining Tn and J with formulae extracted from 
quantum Monte Carlo simulations (see Sec. IIIB, and Ref. 7). The dash in the p3 column for Ag(pyz)2(S208) reflects 
the fact that there is no third component in Eq. (22). Dashes in the J/Zcb and J±/ J columns for [Cu(pyz)2(pyo)2](BF4)2 
indicate a lack of pulsed-field data for this material. 



(e.g. Eq. (3)) is not necessesary in fitting these spec- 
tra, and is set to zero. Fitted parameters are shown 
in Table V. 

Data for the X~ = PFg compound was subject to 
similar analysis, fitting spectra below Tn to Eq. (23), 
this time with pi ^ 15%, ui{T = 0) ^ 12.0 MHz; 
P2 ^ 10%, iy2{T = 0) ^ 9.3 MHz; and ps ^ 75%, 
A3(T = 0) ^ 100 MHz. The phase again proved 
unnecessary. These spectra do not show as sharp a 
transition as the X~ = SbF^ compound, with the os- 
cillating fraction of the signal decaying rather before 
the appearance of spectra whose different character 
indicates clearly that the sample is above T^. The 
available data do not allow reliable extraction of crit- 
ical parameters, but we estimate 5.5K<Tn<6.2K 
and 0.15 < /3 < 0.4. Fitted values are shown in Ta- 
ble V. 

Another method to locate the transition is to ob- 
serve a transition in the amplitude of the muon spectra 
at late times to observe the transition as a function 
of temperature from zero in the unordered state to 
the '|-tair characteristic of LRO, described in Sec. H. 
Spectra were fitted with the simple relaxation function 
A{t > 5 ps) = A\)ge~^^^^, and then the amplitudes ob- 



tained fitted with a Fermi-like step function 

Ait>5^s,T)=A, + -^J^-^^, (24) 

which provides a method of modelling a smooth tran- 
sition between Ai = A{T < Tn) and A2 = A{T > 
Tn). The fitted amplitudes and Fermi function are 
shown in Fig. 19 (c). The fitted mid-point Tmid = 
6.4 ± 0.1 K, and width w = 0.3 ± 0.1 K. Spin re- 
laxation peaks just above Tn, and so one would ex- 
pect that Tn lies at the lower end of this transition. 
Thus, the p+SR analysis suggests Tn = 6.1±0.3K 
(i.e. Tmid — w ±w). This is consistent with the esti- 
mate from v'i{T) and the value Tn = 6.2 K obtained 
from heat capacity 

Members of the M = Ni family exhibit Fp oscilla- 
tions rather like their copper counterparts, with a sim- 
ilar fraction of the muons in sites giving rise to dipole- 
dipole interactions. The results of these fits are shown 
along with those from Cu compounds in Table H. Be- 
cause the nickel data were measured at temperatures 
different from those of the copper compounds and, as 
described in Sec. HID, the muon-fiuorine bond length 
in the compound is sensitive to changes in tempera- 
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Figure 18. Example data and fits for M = Ni magnets. From left to right: (a) and (d) show sample asymmetry spectra 
A{t) for T < Tn along with a fit to Eq. (3); (b) and (e) show frequencies as a function of temperature [no data points 
are shown for the second line because this frequency 1^2 was held in fixed proportion to the first, (see text)]; and (c) 
and (f) show relaxation rate A3 as a function of temperature. In the z^(T) plot, error bars are included on the points but 
in most cases they are smaller than the marker being used. The inset to (f) shows 1^2 = plotted against A3 = A on a 
log-log scale. The black line shows A = z/^ . 



X 


Tn (K) 


ui (MHz) 


1^2 (MHz) 


A3 (MHz) 


Pi P2 P3 P a 


PFe 


6.0(4) 


12.0 


9.3 


90 


15 10 75 0.25(10) 2.8(2) 


SbFe 


12.26(1) 


12.3(1) 


8.98(1) 


80 


25 10 65 0.34(4) 3.1(1) 



Table V. Fitted parameters for M = Ni magnets. Errors shown are statistical uncertainties on fitting and thus represent 
lower bounds. Errors on Ni...PF6 could not be estimated due to the fitting procedure (see text). 



ture, care must be taken when comparing these val- 
ues to those of the Cu family. Linearly interpolating 
the bond lengths for [Cu(HF2)(pyz)2]SbF6 at 9 K and 
29 K to find an approximate value of the bond length 
at 19 K yields f(T = 19 K) = 0.1062 ± 0.0002 nm 
(where the error represents a combination of the sta- 
tistical errors on the fits and variation recorded in 
thermometry, and is thus a lower bound), very similar 
to that measured for [Ni(HF2)(pyz)2]SbF6, r^, f(T = 
19 K) = 0.1068 ±0.0004 nm. 

In spite of being isostructural to the 
[Cu(HF2)(pyz)2]X systems, the dimensionality 
of these Ni variants is ambiguous. Susceptibility 
data fit acceptably to a number of models, and ab 
initio theoretical calculations are suggestive of one- 
dimensional behavior, dominated by the exchange 
along the bifluoride bridges. This is discussed more 
fully in Ref. 16. 



VIII. DISCUSSION 

Fig. 20 collects the results from this paper and 
shows how isolation between two-dimensional layers 



varies over a variety of systems; those presented in 
this paper, molecular materials studied elsewhere, and 
inorganic materials. The primary axis is the experi- 
mental ratio Tn/J. Also shown are the ratios of the 
in-plane and inter-plane exchange interactions, J±/J, 
and the correlation length at the transition, ^(T^j/a, 
extracted from fits to quantum Monte Carlo simu- 
lations. Of the materials in this paper, the least 
anisotropic is [Cu(pyo)6](BF4)2, whose low transition 
temperature is caused by a small exchange constant 
rather than particularly high exchange anisotropy. 
The members of the [Cu(HF2)(pyz)2]X family with 
octahedral anions have k^T^/J ^ 0.33, rather more 
than the k^T^/J ~ 0.25 shown by their counter- 
parts with tetrahedral anions (as has been noted 
previously^). This makes the latter comparable 
to highly 2D molecular systems [Cu(pyz)2] (0104)2 
and CuF2(H20)2(pyz), and the cuprate parent com- 
pound La2Cu04. Below this, [Cu(pyo)2(pyz)2](PF6)2 
exhibits a ratio k^T^/J ^ 0.21(1). The proto- 
typical inorganic 2D system Sr2Cu02Cl2 exhibits 
kBT^/J = 0.177±0.009, and thus l{T^)/a = 280±90. 
The most 2D material investigated in this paper, 
Ag(pyz)2(S208), has k^T^/J = 0.148 ± 0.006, im- 
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Figure 19. The asymmetry at late times A(t > 5 ps) mea- 
sured in Ni(HF2)(pyz)2PF6 along with a fit to Eq. (24). 
Overlaid are regions corresponding to the range of values 
found by fitting the frequency of oscillations as a function 
of temperature with Eq. (4) (5.5 K < Tn < 6.2 K), and 
the value of the A-like anomaly in specific heat capacity 
(SHC) measurements with a representative error of 0.1 K 
(6.1 K < Tn < 6.3 K). The darker area represents the 
overlap between these regions. 



plying ^(TN)/a = 1000 ± 300, with the added benefit 
that the magnetic field required to probe its interac- 
tions is far closer to the range of fields achievable in 
the laboratory. By these measures, molecular magnets 
provide some excellent realizations of the 2DSLQHA, 
with Ag(pyz) 2(8208) being the best realization found 
to date. 

Another method of examining the dimensionality 
of these systems is to consider their behavior in the 
critical region. The critical exponent /3 is a quantity 
frequently extracted in studies of magnetic materials, 
and it is often used to make inferences about the di- 
mensionality of the system under study. In the critical 
region near a magnetic transition, an order parameter 
identical to the (staggered) magnetization, would 
be expected to vary as 



^(T) = $0 1 



T 



/3 



(25) 



In simple, isotropic cases, the value of p depends on 
the dimensionality of the system, (i, and that of the 
order parameter, D. For example, in the 3D Heisen- 
berg model ((i = 3, I) = 3), /3 = 0.367, whilst in 
the 2D Ising model {d = 2, D = I), p = \. Since 
the muon precession frequency is proportional to the 
local field, it is also proportional to the moment on 
the magnetic ions in a crystal, and can be used as an 
effective order parameter. However, Eq. (25) would 
only be expected to hold true in the critical region. 
The extent of the critical region (defined as that re- 
gion where simple mean-field theory does not apply) 
can be parametrized by the Ginzburg temperature Tq, 
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Figure 20. Quantification of the two-dimensionality of the 
materials in this paper and comparison with other no- 
table 2DSLQHA systems. Filled circles show materials 
investigated in this paper; open circles show other molec- 
ular materials; and filled triangles show inorganic systems. 
The low dimensionality is parametrized firstly with the di- 
rectly experimental ratio Tn / J, and then with predictions 
quantum Monte Carlo simulations for both the exchange 
anisotropy, J±/ J [Eq. (6)], and the correlation length of an 
ideal 2D Heisenberg antiferromagnet with the measured J 
at a temperature Tn [Eq. (7)]. The greying-out of the axis 
for J±/J indicates where Eq. (6) is extrapolated beyond 
the range for which it was originally derived^^. Values of 
Tn/ J for the other materials were evaluated from Refs. 50- 
52. 



which is related to the transition temperature Tc by^^ 

2 

d 



(26) 



where d is the dimensionality, ^ is the correlation 
length and AC is the discontinuity in the heat ca- 
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Figure 21. Critical exponent /3, as commonly ex- 
tracted from Eq. (4), plotted against the experimen- 
tal ratio /cbTn/J, indicative of exchange anisotropy. 
Bright filled circles indicate [Cu(HF2)(pyz)2]X materi- 
als examined in Sec. Ill, whilst darker circles indicate 
other materials studied in this work: [Cu(pyz)2(pyo)2]BF4 
(Sec. IV); [Cu(pyo)6](BF4)2 (Sec. V); Ag(pyz)2(S208) 
(Sec. VI). Open circles indicate other molecular materials: 
CuF2(H20)2(pyz) (Ref. 54); [Cu(pyz)2] (€104)2 (Ref. 8). 



pacity. Quantum Monte Carlo simulations suggest^ 
that AC/kB ~ J±/ J- It follows that for d = 3, where 
AC/kB ^ 1 we have \Tg - T,\ /T, ^ (e/a)-^ giv- 
ing rise to a narrow critical region. In two dimen- 
sions, we have \Tg - T^l /T^ ~ / a)'^ {AC / kB)'^ . 
Anisotropic materials with small J±/J only show a 
small heat capacity discontinuity, while ^/a grows ac- 
cording to Eq. (7). This leads to a |Tg - Td /T^ of 
order 1 for our materials, that is, a larger critical re- 
gion for 2D (as compared to 3D) systems. 

The large critical region in these materials allows 
meaningful critical parameters to be extracted from 
muon data. The simplest method of doing so is to 
fit the data to Eq. (4), as we have throughout this 
study; alternatively, critical scaling plots can be used 
(e.g. Ref. 18), which we have performed, finding the 
results are unchanged within error. Since /3 might 
be expected to give an indication of the dimensional- 
ity of the hydrodynamical fiuctuations in these ma- 
terials, a comparison between extracted /3 and ex- 
change anisotropy parametrized by /cb^n/J is shown 
in Fig. 21. Members of the [Cu(HF2)(pyz)2]X family 
show some correlation between the critical exponent 
and the effective dimensionality but overall, the rela- 
tionship is weak. This is probably because /3, which is 
not a Hamiltonian parameter, is not simply a function 
of the dimensionality of the interactions, but probes 
the nature of the critical dynamics (including prop- 
agating and diffusive modes) which could differ sub- 
stantially between systems. 



IX. CONCLUSIONS 

We have presented a systematic study of muon-spin 
relaxation measurements on several families of quasi 
two-dimensional molecular antiferromagnet, compris- 
ing ligands of pyrazine, bifiuoride and pyridine-A/"- 
oxide; and the magnetic metal cations Cu^^, Ag^+ 
and Ni^+. In each case 11+ SR has been shown to 
be sensitive to the transition temperature Tn, which 
is often difficult to unambiguously identify with spe- 
cific heat and magnetic susceptibility measurements. 
We have combined these measurements with predic- 
tions of quantum Monte Carlo calculations to identify 
the extent to which each is a good realization of the 
2DSLQHA model. The critical parameters derived 
from following the temperature evolution of the p+SR 
precession frequencies do not show a strong correla- 
tion with the degree of isolation of the 2D magnetic 
layers. 

The analysis of magnetic ordering in zero applied 
field in terms of inter-layer coupling J± presented here 
does not take into account the effect of single-ion-type 
anisotropy on the magnetic order. This has been sug- 
gested to be important close to Tn in several examples 
of 2D molecular magnet where it causes a crossover 
to Xy-like behavior. In fact, its influence is con- 
flrmed in the nonmonotonic B-T phase diagram seen 
in [Cu(IIF2)(pyz)2]BF4. It is likely that this is one 
factor that determines the ordering temperature of a 
system, although, as shown in Ref. 11, it is a smaller 
effect than the interlayer coupling parametrized by 
J±. The future synthesis of single crystal samples 
of these materials will allow the measurement of the 
single- ion anisotropics for the materials studied here. 

The presence of muon-fluorine dipole-dipole oscil- 
lations allows the determination of some muon sites in 
these materials, although it appears from our results 
that these are not those that lead to magnetic oscil- 
lations. However, the Fp signal has been shown to be 
useful in identifying transitions at temperatures well 
above the magnetic ordering transition, which appear 
to have a structural origin. The fluorine oscillations 
hamper the study of dynamic fluctuations above Tn, 
which often appear as a residual relaxation on top of 
the dominant nuclear relaxation. It may be possible 
in future to use RF radiation to decouple the influ- 
ence of the fluorine from the muon ensemble to allow 
muons to probe the dynamics. 

The muon-spin precession signal, upon which much 
of the analysis presented here is based, is seen most 
strongly in the materials containing Cu^+ and is more 
heavily relaxed in the Ni^+ materials. This is likely 
due to the larger spin value in the Ni-containing ma- 
terials. This is borne out by measurements on pyz- 
based materials containing Mn and Fe ions^^, where 
no oscillations are observed, despite the presence of 
magnetic order shown unambiguously by other tech- 
niques. In the case of Mn-containing materials mag- 
netic order is found with p+SR through a change in 
relative amplitudes of relaxing signals due to a differ- 
erence in the nature of the relaxation on either side of 
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the transition. It is likely, therefore that muon studies 
of molecular magnetic materials containing ions with 
small spin quantum numbers will be most fruitful in 
the future. 

Finally, the temperature dependence of the relax- 
ation rates in these materials has been shown to be 
quite complex, reflecting the variety of muon sites in 
these systems. In favourable cases these data could be 
used to probe critical behavior, such as critical slow- 
ing down, although the unambiguous identification of 
such behavior may be problematic. 

Despite these limitations on the use of p+SR in ex- 
amining molecular magnetic systems of the type stud- 
ied here, it is worth stressing that the technique still 
appears uniquely powerful in providing insights into 
the magnetic behavior of these materials and will cer- 
tainly be useful in the future as a wealth of new sys- 
tems are synthesised and the goal of microscopically 
engineering such materials is approached. 
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